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We propose an efficient dual boson scheme, which extends the DMFT paradigm to collective 
excitations in correlated systems. The theory is fully self-consistent both on the one- and on the 
two-particle level, thus describing the formation of collective modes as well as the renormalization of 
electronic and bosonic spectra on equal footing. The method employs an effective impurity model 
comprising both fermionic and bosonic hybridization functions. Only single- and two-electron 
Green’s functions of the reference problem enter the theory, due to the optimal choice of the self- 
consistency condition for the effective bosonic bath. We show that the theory is naturally described 
by a dual Luttinger-Ward functional and obeys the relevant conservation laws. 


I. INTRODUCTION 

Strongly correlated electron systems remain one of the 
most interesting subjects in modern condensed matter 
physics. It is hard to treat such systems analytically 
due to the large local and, more importantly, nonlocal 
electron-electron correlations, and therefore the devel¬ 
opment of appropriate computational methods is impor¬ 
tant. Nonlocal correlation effects are very important for 

studying charge -ord erinffl^l and Wig ner-MotiPEl tran¬ 
sitions, plasmorP^ and magnorfi^El modes, antiferro¬ 
magnetic fluctuat ion^l^lhSI and other interesting features 
of such systems. These phenomena are realizable in 
adatoms on semiconducting surfaceJ^M^, or in systems of 
cold atomJisHIll and in graphen^i^. 

Dynamical mean-held theory (DMFTjPsEi] bas be¬ 
come the standard approximation for strongly correlated 
fermionic systems. In this theory, all local correlations 
are treated via an auxiliary impurity problem, i.e., the 
electrons on a site are inhuenced by an effective local 
electronic bath formed by the other electrons. This ef¬ 
fective de script ion becomes exact in the limit of inhnite 
dimensioiPSIII] and captures the formation of Hubbard 
band^^mH and the Mott transitiorP^I. In hnite dimen¬ 
sion, DMFT is an approximation that neglects nonlocal 
correlation effects. In particular, the impact of collective 
modes on the impurity problem is neglected. 

There have been many attempts to go beyond DMFT 
and to incorporate these effects. The Dual Fermion (DF) 
approach has been develmed to take nonlocal fermion 
correlations into accounlF^. This approach is exact in 
the two important limits of large and small local interac¬ 
tion. Similar efforts have been made in the DTA^^ and 
IPplapp roaches and in the recently proposed DMF^RG 
methocp328) Nonlocal fermionic correlations are crucial 
for a description of pseudogap formatiorP^, critical expo¬ 


nents for magnetic phase transitioniSHH ^nd formation 
of flat bands near van Hove singularitjPSl. However, these 
methods cannot describe bosonic degrees of freedom and 
their influence on the auxiliary local model. 

Extended dynamical mean-held theory (EDMFTpUlSSI 
was introduced to include bosonic degrees of freedom into 
DMFT, with the main focus on nonlocal density-density 
interactions. In this approach all fermionic and bosonic 
correlations are treated on the effective impurity level, 
and therefore are local. However, it was realized that 
nonlocal corrections beyond EDMFT are necessary for a 
correct description of strongly nonlocal effects, such as 
plasmons. Even though DMFT and EDMFT are similar 
in spirit, the reduction of correlation effects to their local 
part works much better for fermions than for bosons. 

Thus, EDMFT served as a starting point for theories 
that include further spatial correlations. The hrst exam¬ 
ple of such a theory is EDMFT-bGWl^. EDMFT-bGW 
and, more recently, TRILEXP^ diagrammatically treat 
some nonlocal effects beyond EDMFT. However, conser¬ 
vation l aws are not automatically fulhlled in such ap- 
proachepl^. 

In general, studying the collective excitations in 
strongly correlated systems is challenging. Historically, 
two avenues have been explored. One can start from the 
bare Green’s functions, as in the random phase approx¬ 
imation (RPAjPSH, however the bare Green’s functions 
do not contain any information about the spectral weight 
transfer to the Hubbard bands and correlation effects. A 
theory defined in terms of renormalized Green’s func¬ 
tions, on the other hand, contains all the inform ation 
about the Hubbard bands. However, a consistent! 
description of the collective excitations requires not only 
the renormalized Green’s functions but also the renor¬ 
malized vertices and these are numerically very challeng¬ 
ing to handle. 

Recently, the Dual Boson (DB) approactP^ was de- 
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veloped to address these issues. It applies a transfor¬ 
mation to new degrees of freedom that contain the in¬ 
formation about the Hubbard bands and correlation ef¬ 
fects already in their bare Green’s functions. In this way, 
the DB approach for strongly correlated systems fulfills 
conservation law^, necessary for a correct description of 
plasmon^, while remaining computationally tractable. 
The DB method is a diagrammatic extension of EDMFT 
that can be applied to correlated lattice fermion models 
with local and nonlocal interaction. It allows us to in¬ 
clude spatial fermionic and bosonic correlations beyond 
EDMFT by introducing new dual variables. These dual 
variables are introduced via fermionic and bosonic hy¬ 
bridization functions and that act as the effective 
mean-fields acting on an auxiliary single-site impurity. 
Then, the impurity model serves as a starting point for 
a perturbative expansion. By choosing the hybridization 
functions in an optimal way, this perturbative expansion 
can be simplified. In this work, we will study how the 
bosonic hybridization function should be determined 
to optimally treat the feedback of collective modes onto 
the impurity. 

As in DF, the choice of the hybridization function 
has implications for the dual perturbation theory. We 
will show that the correct choice of A^ removes all local 
two-particle processes from the dual perturbation theory. 
Furthermore, we show that the self-consistent determi¬ 
nation of A(j gives the correct results in some important 
limits and that it satisfies the charge conservation law. 
Finally, we study the physical impact of the choice of A^j, 
showing that the nonlocal charge fluctuations make the 
system more insulating. 

The paper is organized as follows: We start by giv¬ 
ing a short description of the Dual Boson formalism 
in section |lTj Then, in section |III| we turn to the self- 
consistency conditions in DB and related methods. In 
section we discuss the dual Luttinger-Ward func¬ 
tional and the generation of self-energy diagrams. The 
impact of the self-consistency on the self-energy dia¬ 
grams is given in section |Vj In section |VI[ we show that 
the self-consistent DB method satisfies charge conser¬ 
vation requirements. Analytical and numerical results 
regarding the self-consistency condition are given in sec¬ 
tion EIIl 


II. DUAL BOSON FORMALISM 

The DI^^ approach to strongly correlated systems re¬ 
lies on a separation of local and nonlocal correlation ef¬ 
fects. In this section we give a general description of the 
approach, and in the next section we will elaborate on the 
role of the self-consistency condition in DB. More details 
on the derivation and the technical implementation can 
be found in Refs. |37| and IH 

We consider the extended Hubbard model in the Mat- 
subara formalism. To be specific, we restrict ourselves in 
this work to the half-filled extended Hubbard model on 



Figure 1. (Color online) Sketch of the DB formalism. The 
original action with parameters U, Vq and tk is replaced by 
an auxiliary impurity problem with fields A„, A„. The expec¬ 
tation values of this impurity model (gr^, and 

enter a dual theory in terms of G and X. 


a square lattice, with the action 
^[c*, c] = - ^ 

jva juj 

+ X + \ X (1) 

(jl)va (jl)u> 

Here, (cjva-) are Grassmann variables corresponding 
to creation (annihilation) of an electron on site j with 
spin (T and fermionic Matsubara frequency v. rij counts 
the number of electrons on site j, and a; is a bosonic 
Matsubara frequency. A normalization by /3 is implied 
in the sum over frequencies. The chemical potential /r is 
chosen in such a way that the average number of electrons 
per site is one (half-filling). U and V are the on-site and 
nearest-neighbor interaction and t is the hopping integral 
between neighboring sites. We use t = 1 as the unit of 
energy. 

The action ([^ is split into a set of single-site impurity 
problems and the remaining part Srem, 

5'imp[c , c] = ^ T p. 

VG 

+ X X! + A,^]n_ij (2) 

S=Y,S%\c\c\ + S,,^. ( 3 ) 

i 

Here, the hybridization functions and A^j were in¬ 
troduced. For the moment, these are arbitrary, and we 
will study the choice of the hybridization functions in this 
work. 

The DB formalism proceeds by introducing new dual 
degrees of freedom /*, /, (j) via a Hubbard-Stratonovich 
transformation of Srem and integrating out the original 
degrees of freedom c* and c. In this way, the original 
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Figure 2 . Interaction vertices in the dual perturbation theory. 


action Q is subdivided into two separate problems: the 
impurity action ([^ and a dual action 


s[rj; </>] = - E Vk.. 

kZ/'(7 

qct) 

+ v[rj, 4 >]- ( 4 ) 

This dual action has bare propagators 

= [9u^ + ~ ^k] ^ — gu, (5) 

- V,]"^ - Xc.- (6) 

Here Sk and Vq are the Fourier transform of tji and V};, 
g and x Eire the impurity Green’s function and charge 
susceptibility respectively 

5,. =-(ci.c*)i^p, (7) 

~ ~ (Pc^P—oj) imp ’ (^) 


where p = riu, — (n) S^i and (.. denotes the impu¬ 
rity average with respect to the action ([^. The single¬ 
frequency susceptibility X = —Ari=i'_ 2 = 2 ' can be obtained 
from the usual expression of the two-particle correlation 
function iFi 22 'i' = (ciC 2 C 2 /C^/) (hereafter, we use com¬ 
bined frequency-spin subscripts, e.g. 1 = i^i,cri). 

The main philosophy of the DB method is that the 
impurity action can be solved numerically exactly, so we 
should try to put as much of the physics as possible into 
the impurity model. The exact solution of the impurity 
model can be fou nd using a continuous-time quantum 
Monte CarlcpiEsl solver. Only weaker correlations re¬ 
main in the dual perturbation theory, which yelds the 
fully renormalized fermionic and bosonic dual propara- 
tors Gk,y = - (/kj./ki.) and Xq^ = - ((^q^(()_q _^). The 
scheme is illustrated graphically in Fig. [l] 

The dual interaction V contains interaction terms of 
arbitrary order in /*, / and (j). All vertices up to three- 
particle level with at least two fermion lines are shown 
in Fig. We use the notation for the vertex with 
n fermion and m boson lines. Usually, the interactions 
are restricted to the two-particle level, where there is 
a fermion-fermion interaction 7 = 7 ^’'^ and a fermion- 


boson interaction A = 7^4 

772.3.4 = (51525354)"^ [ (C1C2C3C4) - 

(cic:) {C2cl) + (cict) (C 2 c:)], (9) 

7 M ;2 = -(5194X2)"^ [ (Cl 04712) - (C1C4) (712)] 

= (5154X2)"^ E [ (ClC 2 -h 3 C 3 C;) - (C1C4) (C2+3C3)]. 

( 10 ) 


These interaction vertices are related 

7i.4;2 = X 2 ^ E [ 7 l.’E. 3 . 4 ff 2-|-353 “ <51.3^2-1-3.4] • (H) 
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In particular, the fermion-boson vertex is non-zero even 
in the non-interacting system ( 7 ^’*^ = 0). A further sum- 
rule for 7^4 js given in Appendix]^ Exact relations for 
the three-particle vertices are given in Appendix]^ We 
will come back to the issue of restricting dual interactions 
to the two-particle level later. 

To calculate the susceptibility Xq^ = — (nn)^^ of the 
original degrees of freedom, one starts by calculating the 
dual polarization (self energy of the dual bosonic prop¬ 
agator) Hqtj = (A°)“^ — X~^. This is done diagram- 
matically using the expansion of the dual action. Then, 
the lattice susceptibility can be obtained from the dual 
polarization as 


4fq'J = 


X^ 4” XuJ HqLj 


-1 


+ A„ — Uq. 


( 12 ) 


Similar equations hold for the relations between dual 
and lattice fermions. In particular, the Green’s function 
is given by 


G 


-1 

ki/ 


9^ + gu Ski/ gv 


-1 

4“ ^ki 


(13) 


where E is the dual fermionic self energy. 


III. SELF-CONSISTENCY CONDITIONS AND 
BOSONIC SUPERLINE 

The dual interaction functional V can be represented 
by an infinite expansion over the vertex functions of the 
impurity problem. Although the series is formally de¬ 
fined for arbitrary hybridization functions Aj, and A(.j, 
the specific choice affects the convergence of particular 
diagrammatic approximations. Physically, and A^.., 
describe the “effective field” generated by electrons on 
the other sites, felt by electrons in the impurity prob¬ 
lem. The freedom to choose them allows us include most 
of the correlation effects into the impurity, simplifying 
the dual action. To do this, the hybridization functions 
are determined self-consistently according to some self- 
consistency condition. There are several options for the 
self-consistency condition, and here we will discuss them. 
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DMFT is restricted to systems with Vq = 0 and it does 
not use a retarded interaction, = 0. It assumes that 
the self-energy is local, Ski/ = 0- The self-consistency 
condition on A „ is 


According to Eq. (19), the bosonic “superline” is ex¬ 


EGki/ = 


ffi/ 


or, equivalentljf^ in the DMFT approximation, 

5]Gki/ = 0. 


(14) 


(15) 


In both expressions, denotes a momentum average. 

The DF extension of DMFT shines light on the mean¬ 
ing of the self-consistency conditiorP^. The self-energy 
in DMFT is the zeroth-order of the dual perturbation 
theory, and the condition (15) ensures that the first-order 
(Hartree) contribution to the self-energy is zeroIn this 
way, DMFT and DF only differ in higher orders of the 
dual perturbation theory. In DF, the condition (14) is 


not equivalent to (15) due to nonlocal parts of the self¬ 
energy. 


actly the difference between the lattice and impurity 
susceptibilities. It contains both types of two-particle 
processes, the fermionic ladder and the bosonic fluctua¬ 
tions, and combinations of these two. Using Eq. 
the “lattice” self-consistency condition 0 is exactly 
the requirement that the local part of the “superline” 
is zero (see Eq. (20l). On the other hand, the “dual” 


self-consistency condition (181 can be understood as the 


requirement that the local part of only dual bosonic fluc¬ 
tuations is equal to zero. The former self-consistency 
condition treats collective fermionic and bosonic fluctua¬ 
tions on equal footing. 

Until now, moslP^ numerical results obtained with 
Dlj ^ | 8 | 9 l started with EDMFT (Hq^ = 0) self-consistency 
and then did only a single-shot of DB, so the question of 
self-consistency was not relevant. These results showed 
that the susceptibility in DB satisfies charge conserva¬ 
tion. This is a useful property, crucial for studying long- 
wavelength collective excitations, and we make sure that 
our self-consistent calculations will also satisfy this re¬ 
quirement. This issue will be studied in Sec. m 


^ _ ^imp 


1 4” 


(16) 


and Eq. (14) does not make the first-order self-energy di¬ 


agram vanish. This is a strong motivation to use Eq. (15) 
in DF. 

EDMFT, on the other hand, adds a retarded interac¬ 
tion to account for impurity screening by the nonlocal 
interaction Vq. This retarded interaction is determined 
by a self-consistency condition similar to that for 
namely 


A^qu Xo. 


q 


(17) 


Similar to the fermionic hybridization, within EDMFT 
there is an equivalent self-consistency condition (se^ or 
Appendix 


5]Aq^ =0. 


(18) 


Finally, the DB approach also takes into account the non¬ 
local Ilqcj and a choice between (18) and 0 needs to 
be made. 

To study the meaning of the self-consistency condition 
( [T7| ), we use an exact expression for the lattice suscepti¬ 
bility in terms of dual variables (see Appendix [d| 


A = (i + xn)[x(i + nx) + x] =x + ^, 


(19) 


where we define the bosonic “superline” and propose the 
new self-consistency condition as 


q 


^qi. = 0, 


5 = A + xiix + xiiA + xfix + xnAny. (20) 


IV. DUAL LUTTINGER-WARD FUNCTIONAL 

Let us now discuss a functional description of DB, 
containing both collective fermionic and bosonic fluctu¬ 
ations. Before present ing a Baym-Kadanoff functional 
derivation of the theorjB^Sl^ need to analyze how the 
conserving properties of the original and dual theory are 
related. We will show that the origmal and dual theo¬ 
ries are simultaneously <l>-derivable^, so that the self¬ 
energy function fulHll variational equations <5$ = E5G 
and 5$ = Ei5G with certain functionals <i>, $. This is 
equivalent to the statement that 5$ — 5$ is a full differ¬ 
ential. 


Using the notation i? = 1 -|-gS, the exact relations (13) 
and (16) can be expressed in a compact form. 


G = Rg + R^G, 


( 21 ) 


which is the equivalent of Eq. (20) for fermions, and 


T, = w — A — (Rg) 


( 22 ) 


An important property of this representation is that the 
dispersion law does not explicitly enter in the formu¬ 
las. Therefore only renormalized quantities appear in the 
variational procedure, in the spirit of the Baym philoso¬ 
phy. We obtain 

= {ii,-A)6G-{Rg)-^6G-g-\R-l)SG. (23) 

We consider the variational procedure at fixed hybridiza¬ 
tion, so that variances of A and g do not appear. Con- 
cequently, using the relation (21), the calculation gives 


= {w-A)6G-{R-^+2g-^G)6R-g-\2R-l)SG, 

(24) 
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and we can get 

$ - I. = (ii/ - A)G + g-^G - Ini? - 2g-^RG. (25) 

This result is very important, because it shows the gen¬ 
eral possibility to construct a conserving theory from the 
dual theory. The obtained expression also establish an 
explicit relation between the original and dual function¬ 
als. 

Now we can introduce the dual functional for our the¬ 
ory to show how to generate self-energy and polarization 
diagrams. The dual functional formulation does not lead 
to a straig htforw ard proof of the Baym-Kadanoff conser¬ 
vation law^^ 2 llll fQj. original lattice problem. Charge 
conservation will be discussed in Sec. I VII 

7 \ 7 \ = + 


Figure 3. (Color online) Diagrammatic representation of 
Eq. (20l. The dual superline S describes two-particle fluc¬ 
tuations that have both pure bosonic and collective fermionic 
character. Filled lines denote the dual propagators, dashed 
lines the impurity ones. 


Figure]^ is a diagrammatic representation of the “su¬ 
perline” in Eq. (20), where the two-particle ladder was 
inserted for 11 since we only want to consider two-particle 
processes. The l add er is necessary for charge conser- 
vatiorP (see Sec. VI). These diagrams contain only a 
dual bosonic line and two-fermionic ladders, so they de¬ 
scribe all two-particle processes and no processes involv¬ 
ing more than two particles. 



Figure 4. (Color online) Dual Luttinger-Ward functional for 
the Dual Boson approach. The generation of diagrams for the 
dual self-energy and polarization function occurs by cutting 
one fermionic or bosonic line respectively. 

As usual, the diagrams for the fermionic self energy and 
polarization operator can be found diagrammatically by 
cutting the respective line of the associated dual func¬ 
tional (see Fig. 1^, 


n — Q + I ; 

□ = □ + O’ ^ 


< + <□ 



Figure 5. (Color online) The renormalized three- and four- 
point vertex functions in the Dual Boson approach. The five- 
and six-point vertices are introduced here with the local dual 
fermion line. Such diagrams are identicaly equal to zero be¬ 
cause of the self-consistency condition and inserted only to 
generate the self-energy diagrams from the Luttinger-Ward 
functional. 



= <ct> + <m> 


Figure 6. (Color online) Diagrams for the dual self-energy 
and polarization function. The renormalized three- and four- 
point vertices have the ladder structure (see Fig.[^ and don’t 
contain any diagrams with the local fermionic line due to the 
self-consistency condition. 


self-consistency condition (15), there is no need to ac¬ 
count for closed fermionic loops in diagrams for the self¬ 
energy and polarization operator. However, this is not 
the case for the dual functional diagrams. Indeed, cut¬ 
ting a closed loop in the dual functional yields a non-zero 
contribution into the local part of the self-energy. To ac¬ 
count for this, we introduce the renormalized vertices of 
the dual functional. Each renormalized vertex is the sum 
of a bare part and a higher-order vertex with a closed 
fermionic loop (see Fig. [^. These renormalized vertices 
are used as building blocks for the ladders in the func¬ 
tional of Fig. 1^ This is the minimal possible functional 
that describes the full physics of two-particle fluctua¬ 
tions. 

The diagrams for the self-energy and polarization can 
be determined from the functional. These are shown in 
Fig.i with renormalized three- and four-point ladder 
vertex functions as in Fig. 


V. SELF-CONSISTENCY CONDITION AND 
HIGHER-ORDER VERTICES 


11 =^. 
SG SX 


(26) 


There is however an important peculiarity. Since the lo¬ 
cal part of the fermionic propagator is zero due to the 


The dual action contains vertex functions of higher or¬ 
der than the two-particle correlators 7 ^’^ and 7 "^’° (see 
Fig. §. There is no small parameter in the dual per¬ 
turbation theory, so in principle the contribution of all 
diagrams containing these additional vertices should be 
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taken into account. However, this is usually not feasible 
in practice. Neglecting higher- order vertices can lead to 
small but noticeable deviation^^SMl, Although they are 
usually small, it would be good to reduce these deviations 
further. 

The fermionic self-consistency condition already re¬ 
moves all diagrams with higher order vertices that con¬ 
tain local dual Green’s functions. In the dual fermion 
approach, the hybridization function Aj^ is the only free 
parameter, so there is only one self-consistency condi¬ 
tion, and there is no more freedom to simplify the dia¬ 
grammatic expansion. However, the DB approach has an 
additional free parameter, the retarded interaction 
and it can be chosen self-consistently in such way as to 
minimize the impact of higher-order vertices. 

We recall that DMFT captures the single-electron hy¬ 
bridization physics and is formulated in terms of the 
single-particle Green’s function. The DMFT choice of 
Ai, removes the lowest-order Hartree diagrair^ contain¬ 
ing the two-particle quantity 7 ^’°. Here we would like to 
include the bosonic (or, equally, two-fermionic) effects. 
In analogy to the foregoing, we construct a theory in 
terms of single- and two-particle quantities, and require 
that the lowest-order diagrams with the three-particle 
quantities 7 ®’°, 7 ^’^ and 7 ^’^ drop out of the self-energy. 
These lowest-order diagrams are shown in Fig. where 
E”’™ denotes the contribution to the dual self energy 
from the 7 "’”^ vertex function. 

These contributions to the dual self-energy from 
higher-order vertices can be divided into reducible and 
irreducible part^^. Here we consider the irreducible con¬ 
tribution from the higher-order vertices; Appendixj^de- 
scribes the reducible contributions. 

The three-particle vertex 7 ®’° in these diagrams is 
connected to four internal fermion lines. Gounting the 
number of adjustable parameters, it is easy to realize 
that a four-time (three frequency) bosonic hybridization, 
Ai 22 ' 1 'c|c 2 C 2 'Ci' is needed to completely remove the dia¬ 
grams of Fig. 

Such a quantity is very hard to implement in any prac¬ 
tical calculation. Therefore, we stay with the two-time 
(one frequency) hybridization A„. In this way, we can¬ 
not make the contribution from the diagrams with three- 
particle correlators exactly zero. However, we will show 
that the self-consistency condition ( [^ corresponds to 
removing the contribution from the physically important 
part of 7 ®’®. In the case of small [/, all three-particle cor¬ 
relation effects vanish because of the perturbative form of 
the diagrammatic expansion of the dual functional (see 
Appendix . 

The three-particle vertices in the dual perturbation 
theory are impurity correlation functions. Since the im¬ 
purity only contains two-particle interaction [see Eq. ([^], 
7 ®’® can be written as a set of diagrams containing 
the irreducible two-particle vertex function of the impu¬ 
rity problem connected via local fermionic and bosonic 
Greens functions. Such a diagrammatic expansion of the 
dual-fermion six-leg vertex was discussed in Ref. 1501 



Figure 7. (Color online) The structure of the irreducible parts 
of higher-order vertices with respect to the bosonic argu¬ 
ments. The diagrams a) and b) give the contribution S®’® 
to the self-energy from the six-point vertex, diagrams c) and 
d) are the contributions and E^’'^ from the vertices 
and 7 "^’^ respectively. The purple parts of the diagrams are 
equal to the dual polarization function H (see Fig. 1 ^, the 
dashed wave lines are the impurity susceptibility y. 


Here we analyze three-particle dual self-energy dia¬ 
grams where the incoming particle is connected to a two- 
particle ladder, which depends on a single bosonic fre¬ 
quency. Fig. shows these processes. Physically, these 
diagrams describe the interaction of a fermion with a col¬ 
lective excitation. The collective excitation can take the 
shape of a dual boson propagator or of a ladder of dual 
fermion propagators. 

Let us start with diagrams a) and b). These diagrams 
give the following contribution to the self-energy, where 
we denote the purple parts of the diagrams in Fig. as 
the dual polarization function H: 

E“ + S^ = 172.2^ [n +HAH] y. ( 27 ) 


Here | appears as a symmetry factor due to indistin¬ 
guishable vertices. 

The contributions of diagrams c) and d) can be written 
in the same way. 


E" = 

E'^ = (28) 


where | in S'” is also a symmetry coefficient. 

Then, the total contribution to the dual self-energy 
from these diagrams is 


E® = E“ 


■ E^ -b E® -b E'^ 


=- 72^2 + 2Anx + x[n +HAH] y). (29) 


Here, we recognize the superline of Eq. (20). This allows 
to rewrite (29) as 




(30) 


where the local part appears since both ends of the su¬ 
perline are attached to the same vertex. This expression 
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is equal to zero if we use the self-consistency condition 
(20). Therefore, the local part of the “superline” S being 


zero means that the dual self-energy contribution asso¬ 
ciated with 7 ^’^ drops out of the perturbation theory. 


For comparison, the self-consistency condition (18) 


corresponds to = 0. In that case, diagrams a), b) 
and d) still contribute. 

Having removed the lowest-order three-particle pro¬ 
cesses with three-particle vertices, we can now restrict 
ourselves to the the simple ladder expansion for the dual 
self-energy and polarization function 


S = 


n = . 


+ 







(31) 

(32) 


The reducible contributions to the three-particle vertexes 
can be accounted by the re-definition on the lines in these 
diagrams, as it is described in Appendix]^ These expres¬ 
sions allow to close the set of equations for dual Green’s 
functions formed by the self-consistency conditions (15 
20 ) and definitions ([^ . 

Thus, we obtained a diagrammatic argument for the 
bosonic self-consistency similar to that in the DF ap¬ 
proach. This condition removes the physically important 
dual self-energy contribution from higher-order vertices 
and simplihes the diagrammatic expansion of the inter¬ 
action functional V. The main difference is that for the 
bosonic sector, both single-boson and two-fermion con¬ 
tributions are crucial. 

It is important to note that a nonlocal self-energy can 
break the charge conservation law, since the vertex ob¬ 
tained as a second derivative of the dual Luttinger-Ward 
functional then depends on three momenta, whereas the 
ladder approximation only takes one transferred momen¬ 
tum into account. For this reason, in this work we 
restrict ourselves to a local self-energy and focus on 
nonlocal bosonic effects of Eq. (32). As a result, the 


two self-consistency conditions (14) and (15) for the 


fermionic hybridization are equivalent. The fermionic 
and bosonic self-consistency conditions remove the main 
physical contributions to E and the fermion physics is 
accounted for on the level of the impurity. 

Our computational scheme is summarized in Figure]^ 
We start with an initial guess for the hybridization func¬ 
tions and A„. In this work we take the result of a con¬ 
verged EDMFT simulation as the inital guess, however 
it is also possible to start from non-interacting Green’s 
functions or from a previously obtained converged DB re¬ 
sult at other parameters. The impurity problem is solved 


using quantum Monte GarlcPSHSEIMl and the resulting 
impurity quantities are used to calculate the dual polar¬ 
ization. Then, the lattice Green’s function and suscepti- 
bilty are determined and new hybridization functions are 
determined using the self-consistency conditions. This 
procedure is iterated until convergence is achieved. 


Ai/, A^^ 


O 

CO 


Gk., A, 


qo; 


impurity solver 

- > 

''' 

ly','' 

if' 


( 12 ) 


^ 2.1 




Figure 8. A summary of the computational scheme used. The 
impurity solver determines local impurity quantites based on 
hybridization functions A^,A^. The fermionic lattice Green’s 
function is determined directly from the impurity quantities. 
For the lattice susceptibility, we take additional nonlocal cor¬ 
rections into account via the dual polarization ll in the ladder 
approximation. 


VI. CHARGE CONSERVATION 


Here, we prove that the susceptibility obtained from 
DB satisfies the charge conservation requirement as long 
as the bare (local) self-energy and the local vertices are 
used. The conditions for the proof are: 

• The impurity problem ([^ is solved exactly, so that 
the impurity vertex and self-energies satisfy the 
Ward identity. 


• The self-energy of the fermions is local and equal 
to the impurity self-energy, Sk,y = E™p. 

• The hybridization is chosen in such a way that 
the local part of the dual fermion Green’s function 
is zero, i.e., the self-consistency condition (15) is 
satisfied. 


It is important to note that for this proof there is no 
condition on the retarded interaction A^^. In particular, 
the proof works for the self-consistent choice of A;^ in the 
previous section. 

We start with the dual polarization in the ladder ap- 
proaclP, Eq. (32), and write it as 






p2,l 


(33) 


where we have defined the bubble of dual Green’s func¬ 
tions 


Gk!^Gk-t-qZ/-t-Cd5 


(34) 


and renormalized fermion-boson and fermion-fermion 
vertices 


r2.i = ^2,1 _ py.o ^ 

qi/'o; luu) /j aui''uj ' 


ciiy'uj !jy>^ 


y 4,0 _ 4,0 _ 

' UJ ~ 2 ^^ > 


4.0 


.B. 


^4,0 


uv" c]y"v'uj’ 


(35) 

(36) 


























The results of the ladder summation can be expressed as 
a geometric series, 


p 2 ,l _ 


r 


i:[ 

v" 


I+ B, 


4,0 




-1 


-1 


2,1 


4,0 

VijJ ' 


(37) 


where the inverse should be understood as a matrix in¬ 
version in frequency space and / is the identity matrix in 
this space {w"). 

Now, we proceed by going from the impurity vertex to 
the particle-hole irreducible impurity vertex (see Refs. E^ 

and [S3]). 


E r J I Riiiipp4,0 ir 

ViV2V3Vi 


I -1 


c{yz(jj 


-^4,0 ir 






-1 


— Bc^viuj [i + r 


i4,0ir p 


-1 


( 41 ) 


In the Hubbard model {V = 0), this is the expression for 
the DMFT susceptibility. Here we have introduced the 
lattice bubble 


Bq... =BZ^ + Bc^uu. 


J^G'k.Gk+q, 




(42) 


BZ^ = - 9i^au+u^, 


4,0 




1 -1 


Xc 




.4,0 ir 

-1 


(38) 

(39) 


The impurity susceptibility Xuj is related to this irre¬ 
ducible vertex as 




n -1 


^imp 


(40) 


Now, all relevant two-particle correlators of the im¬ 
purity have been expressed in terms of two quantities, 
the irreducible vertex and the particle-hole bub¬ 

ble H™P. This allows us to write the combination 
P = —Xuj~XuiBqi^Xui that enters the lattice susceptibility 
in Eq. (12) entirely in terms of these two quantities. 


This relation holds, since Gki/ = <?;/ + Gki/ and the local 
part of the dual Green’s function is zero so cross-terms 
gG vanish in the momentum sum. 

For the charge vertex 


p4,0 ch 




(43) 


we have the Ward identitj^^ 


= Gk,t, - G^E. = -(^ - S.+. + S.). (44) 


Therefore, for w 7 ^ 0, 


-.4,0 ch 


Pq=0uj — E! q=0l/t.J 

1 / 

= ~ E E Gki'Gk,j/-i-i^[Gj^ ,2 ~ P’Vi.v+uj] 

V k 

= E! E![Gkiy ~ Gk.jy-l-t,;] 


^qu; — Xg; Xgj^ciujXuj 

= -V - Tj I Dimpp 4 , 0 ir 

Xoj / ^ ^ 1/2 ^" 2 ^ 


Bn 


-1 


I + V + ■ 


- -X^- E [l + Bl^SB 


4,0ir 
1 / 21^1 u;. 


X [i + rtfZZBZS + Bq.,.) 

-1 


1 / 11 / 40 ; V 1 / 40 ; 

[r I Rimpp4,0i 

/ > L ~ 1^21^1 

Z/ 11/2 2^4 

X [bzSS.... 

rr, oimpp4,0i 

/ > L “ 2/2 2^2 2^1 


J40 


-1 


-1 






1 -1 


1 / 12 / 40 ; q2^4^j 

-1 


BZ^ [I + TlZZBq..^ + Bq.,.^. 


X [/ + nzi.B, 


4,0 ir p 1 ^ 

2/3 


= 0 . 

The physical polarization is given by 


TT”^ — P~^ — A 

^^qo; ~ ^qcj 


(45) 


(46) 


which yields the required condition nq=o,i.j>o = 0. For 
arbitrary q the physical susceptibility is given by 


/fqa; — 


1 


1 


-Hq^i -Vq-U 


Pqu) -|- G -|- Hq — 


(47) 


which gives 


X, 


q= 0 ,a; 7^0 


= 0. 


(48) 


So the charge conservation requirement is satisfied. 


VII. APPLICATION OF THE 
SELF-CONSISTENCY CONDITION 

Having introduced a new self-consistency condition 
( |17[ ), we will study its behavior in some important lim¬ 
its. In the limit of infinite dimensions, DB with the new 
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self-consistency condition reduces to the (exact in d = 
oo) DMFT. In the non-interacting limit, the new self- 
consistency condition is automatically satisfied by the 
exact solution. At weak interaction, the self-consistency 
condition changes the effective impurity interaction. Fi¬ 
nally, we determine the high-frequency limit of the re¬ 
tarded interaction and find that it is in general nonzero. 


to the self-consistency condition on the fermions, and as 
a result the local part of the bubble [GGJq.o; is also zero. 
So the local part of the lattice susceptibility (fl^ is 


(51) 

q q q 

= X-, (52) 


A. Limit of infinite dimensions 


and the self-consistency condition 0 is satisfied. 


As is well known, DMF T becomes exact in the limit 
of infinite dimension^^HIlI] Qur DB scheme inherits this 
property. To prove this, it is only necessary to show 
that DB and DMFT are equivalent, i.e., that the self- 
consistency condition 0 is satished for vanishing re¬ 
tarded interaction in the limit of infinite dimensions. 
Since nonlocal interactions only contribute at the Hartree 
level in infinite dimension^^, it is sufficient to prove this 
for t/q = 0 . 

w I ft 

Aqc^; — \uj T Xt.iJ-l-lqiA;X£A; 

^ ' Alq,^ ^ ( ^quiXi^ 

q q 

q q 

= 0, (49) 

since the local part of the second-order diagram van¬ 
ishes due to the fermionic self-consistency (as before), 
and the local part of higher-order vertex corrections 
vanishes in infinite dimension^^. The vanishing of the 
higher -order diagram can also be understood in a 1 /d ex- 
pansioiPSED; the four dual Green’s functions contribute 
0{l/\/d) eaclP^, and there is a single internal summation 
contributing 0{d), so llq^^^ oc 0(1/d). 

Eq. (49) shows that the self-consistency condition 0 
is automatically satisfied for A,^ = 0 in infinite dimen¬ 
sions, so DB and DMFT are equivalent (and both exact) 
in this limit. 


C. Strong and weak interaction 

In the special case of the Hubbard model (that is, 
V = 0), the interaction is purely local, but the collective 
modes are still present. Magnons are known to determine 
the low-energy physics of the half-filled Hubbard model. 

The dual fermion ladder summation allows us to de¬ 
scribe the effect of magnon fluctuations, leading to the 
formation of an “antiferromagnetic” pseudogap slightly 
above the Neel temperature^. In the same vein, one 
should expect the presence of collective density modes. 
Indeed, it is widely accepted that the sligtly doped Hub¬ 
bard model shows instabilities related to the formation 
of static density waves. This indicates the importance of 
dynamical density fluctuations in a wider region, possibly 
including the case of half-filling. 

In the dual-boson theory, the presence of collective 
modes should be reflected in a corresponding bosonic hy¬ 
bridization A. Indeed, our self-consistency conditions 
impose that A is non-zero even for V = 0. In this case 
the bare bosonic line is purely local, and the nonlocality is 
provided by the fermionic ladder. The latter exactly cor¬ 
responds to a collective mode, as it was discussed above. 

Two important limits of the Hubbard model (V = 0) 
are those of small and large U. Both of these limits are 
described by Xt^^w “C I: for [/ = 0 the A^^ = 0 as it 
was shown in Sec. |VHBl for [/ — >■ 0 the A^ —0 while Xuj 
stays finite. Conversely, for large U, the susceptibility Xuj 
becomes very small. This allows us to expand the charge 
susceptibility ( 12 ) in Xuj^ui 


B. Noninteracting system 


Alqo; ~X^ 4” X(-^j X(^^(^Xi^ 4- ... (53) 


The noninteracting system (U = 0, V = 0) is ex¬ 
actly solvable, the bare Green’s function is the exact lat¬ 
tice Green’s function. This result is recovered in self- 
consistent DB if A^ = 0. Here, we will show that the 
new self-consistency condition is indeed fulfilled when 
U = 0 and A^ = 0. Starting with a noninteracting 
impurity problem gives = 0 , and consequentlj®^ 

= —1/xuj- Only the second-order diagram 

x^fi^qlx^ = [GG]q,^ (50) 


The polarization also simplifies in the limits of both 
strong and weak interaction. With the same argument 
as in DfUll, the vertex 7 is small at small U, whereas 
the dual Green’s function G is small at large U. This 
means that both at large and at small U, higher orders 
in the ladder decay quickly, and when combining the self- 
consistency condition 0 and Eq. ( [^ , only the third- 
order contribution llO) needs to be considered (the local 
part of the second-order diagram drops out due to the 
self-consistency condition on the fermions, as before). 


contributes to the dual polarization, with G the Green’s 
function of the dual fermions. Its local part vanishes due 


Xu^^u^Xu, ^'^XiJ^qlxuj- (54) 

q 
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D. Symmetries in the Hubbard model 


At small [7, the third-order polarization diagram is 
proportional to U, since it contains a single vertex 7. In 
this work, we focus on the charge sector, however, there 
are also collective excitations in the magnetic sector and 
the DB formalism can be also applied to SzSz interac¬ 
tions. An important difference between the charge and 
magnetic channel is the sign of the effective interaction. 
This implies, to lowest order in U, opposing signs in the 
third-order dual polarization 

and effective impurity interaction -j-AJJ^^snetic _ 

By using particle-hole symmetry, one can show (see, 
e.g., chapter 11 of Ref. 1551 for a pedagogical discussion) 
that the transformation U ^ —U in the half-filled square 
lattice Hubbard model interchanges charge and spin fluc¬ 
tuations. Single-particle properties like the self-energy 
are invariant under the the transformation U —>■ —U. In¬ 
troducing a self-consistency only on the charge sector, 
as we do here, will break this symmetry. To retain the 
?/—>■—[/ symmetry, the charge and magnetic channel 
need to be treated on the same footing. To also retain 
spin rotation symmetry, and in fact the full sym- 

metrjiSH^ all three magnetic channels are needed, as well 
as the charge channel and two superconducting channels. 


E. Screening by nonlocal interaction 

Nonlocal charge fluctuations play an especially impor¬ 
tant role in systems with nonlocal interactions. The 
nonlocal interaction can cause quantitative differences by 
screening the effective local interaction. However, it can 
also lead to qualitatively different physics. The checker¬ 
board order that arises from repulsive nearest-neighbor 
interaction is a good exa mple of this, that has been stud¬ 
ied both in EDMFTj^lSZI and in dual bosorP. Attractive 
nearest-neighbor interaction, on the other hand, can lead 
to phase separation into high and low density phase^^. 
Clearly, the sign of the nearest-neighbor interaction V 
is physically relevant. Using the Bogoliubov inequality, 
screening effects have been estimated to be proportional 
to Vl^. 

A peculiarity occurs for EDMFT on the square lattice. 
As shown in Appendix due to the structure of the 
momentum sums, only the absolute value of V matters 
to EDMFT. In dual boson, on the other hand, the sign 
of V does matter. 


F. Numerical results at weak interaction 

To start our discussion of the numerical results, we look 
at a system where correlation effects are expected to be 
moderate. We take a weak local interaction U = 0.5 and 
relatively high temperature T = 1, and vary the nearest- 
neighbor interaction between V = —0.03 and V = 0.03. 
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Figure 9. (Color online) Left: Effective impurity interaction 
as a function of Matsubara frequency for various nonlocal 
interaction strengths V, a.tU — 0.5. The dashed line indicates 
the high-frequency asymptotic value of A,^, and this value is 
shown in the top right as a function of V (the line is a linear 
fit for \V\ < 0.01.). Center-right is the zero-frequency A,,, 
and bottom right the zero-frequency A„ in EDMFT. 


All parameters are given in units of t = 1. The effec¬ 
tive impurity interaction after the dual self-consistency 
is given in Fig. 

One thing that is immediately noticeable from Fig. 
is that the high-frequency asymptotic value of the ef¬ 
fective impurity interaction depends on the nonlocal in¬ 
teraction strength. This is unlike EDMFT, where the 
asymptotic is given by the local interaction strength only 
(see e.g. Fig. 9 of|33or Fig. 6c of [ST]). Fig.[^ shows 
that the asymptotic value depends linearly on V. This is 
reminiscent of the effective local interaction U* = U — V 
obtained using the Peierls-Feynman-Bogoliubov varia¬ 
tional principl^^ where V is proportional to V. The 
linear fit in Fig.[9jD gives V « 1.9U in this case. 

Alternatively, the zero-frequency part of the retarded 
interaction can be used to estimate the effective inter¬ 
action. As mentioned before, in EDMFT this leads to 
a quadratic dependence of the effective interaction on 
V. We have observed the same quadratic behavior in 
EDMFT, see Fig. Hji. The zero-frequency part of the 
retarded interaction in self-consistent DB is shown in 
Fig. [^. We find that the zero-frequency part of the 
retarded interaction also depends linearly on V, with 
V « 0.3U. This constant of proportionality is almost 
an order of magnitude lower than that of the high- 
frequency retarded interaction. On the other hand, the 
zero-frequency retarded interaction in DB (Fig.Om is an 
order of magnitude larger than in EDMFT (Fig. |M). 

Fig.shows that even at U = 0 (Hubbard model), the 
self-consistency condition introduces a retarded interac¬ 
tion to the impurity problem. In the absence of a non¬ 
local interaction, the susceptibility still has a nontrivial 
momentum-dependence and the impurity feels the effect 
of these nonlocal charge fluctuations. An important re¬ 
mark here is that, according to Sec. |VH C[ the retarded 
interactions in the charge and magnetic channel cancel 
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Figure 10. (Color online) Left: Effective impurity interaction 
as a function of Matsubara frequency for various nonlocal 
interaction strengths V, at [/ = 10. The dashed line indicates 
the high-frequency asymptotic value of A^j, and this value is 
shown in the top right as a function of V (the line is a linear 
fit for \ V\ < 0.1.). Center-right is the zero-frequency and 
bottom right the zero-frequency A^i in EDMFT. 


to first order in U in the Hubbard model. So the finite 
retarded charge interaction is mostly compensated by the 
magnetic interaction. 

As an aside, the charge conservation requirement 
-Aq=o,tcj>o = 0 is satisfied numerically in these self- 
consistent DB calculations (not shown). 


G. Numerical results at strong interaction 


Now, we turn to a strongly interacting system, with 
U = 10, and local interaction varying between V = —0.3 
and V = 0.3. The retarded interaction is shown in 
Fig. As before, there is a constant retarded interac¬ 
tion in the hig h-fr equency limit, the linear fit of this co¬ 
efficient in Fig. 10] gives V « 1.9V. The zero-frequency 
retarded interaction also depends linearly on V with co¬ 
efficient V « 0.2V. 

Note that the retarded interaction is quite large and 
positive (repulsive), so the effect of the nonlocal charge 
fluctuations is to make the system more insulating. 

The retarded interaction has a peculiar shape, with a 
strong repulsion at finite frequency and a much smaller 
value at zero frequency. The retarded interaction is based 
on the lattice susceptibility, which has a discontinuity at 
q = 0 and a; = 0. In that sense, it is not completely sur¬ 
prising that the retarded interaction also changes sharply 
around zero frequency. However, one should keep in mind 
that the results shown here only consider fluctuations in 
the density channel. As explained in Sec. VH D[ a full 
computation also needs to take the magnetic and super¬ 
conducting fluctuations into account and we expect these 
fluctuations to be very important for strong interaction 
strengths. 

The change in the impurity interaction affects the lo¬ 
cal observables obtained from the impurity. As an exam- 


Figure 11. (Color online) Quasiparticle renormalization fac¬ 
tor Z as a function of the nonlocal interaction strength V, at 
U = 10. Nearest-neighbour repulsion (V > 0) screens the lo¬ 
cal interaction and increases the quasiparticle renormalization 
factor. 


pie, in Fig. [g we show the quasiparticle renormalization 
factor Z = {1 — Im Eyj^/izi)”^ determined from the local 
self-energy at the first Matsubara frequency. The lattice 
self-consistency condition takes into account feedback of 
the nonlocal charge correlation effects onto the impurity. 
These correlations drive the system in the direction of the 
insulating phase and this is reflected in the reduced Z- 
factor. Secondly, a repulsive nonlocal interaction V > 0 
screens the local impurity problem, makes the system less 
correlated and increases the quasiparticle weight. Within 
the EDMFT approximation, the dependence on E is a lot 
weaker. 

The charge conservation is also visible numerically. 
The susceptibility in single-shot DB at the end of 
EDMFT self-consistency is known to satisfy the charge 
conservation requiremeniP A'q=o,w>o = 0- This is shown 
in Fig. as iteration 0. Now, when the additional DB 
self-consistency is started, and A^^ is updated, this charge 
conservation is initially broken, as can be seen in itera¬ 
tion 1 and 2 of Fig. [I^ Charge conservation is restored 
as self-consistency is achieved. 


H. High frequency asymptotic of 

In the previous section, it was shown that the retarded 
interaction A^^ goes to a constant at high frequency. This 
high frequency asymptotic behavior can be understood 
in the following way. We start by examining the ingredi¬ 
ents of the self-consistency formula. The impurity sus¬ 
ceptibility Xcj decays as At high frequency, the 

fermion-boson vertex increases as cx oc A 
bubble of dual Green’s functions decays as so 
the second-order diagram Hq^L oc XuiBq^ui^ui cx uP'. This 

results in a constant x^^nq^i,, which is necessary for the 

~ (' 2 ) 

charge conservation condition XwnqAQ^_,,gQ = — 1. 

So the magnitude of the second-order diagram is in- 
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Figure 12. (Color online) Charge conservation requirement 
^q=o,t^>o = 0 during the self-consistency scheme. Iteration 
0 is single-shot DB. These results are a,t U = 10, V = 0.3. 


creases as a function of frequency. However, here we only 
consider its local part, which vanishes since the local part 
of the bubble of dual Green’s functions is exactly zero due 
to the fermion self-consistency condition, i.e., 

En5 = 0- (55) 

q 


The higher order ladder diagrams have additional rungs 
Since 7 goes to a constant ({ 7 ™purity^ large 
frequency, every additional order in the ladder decays by 
an extra factor coming from the bubble. That means 
that the first non-vanishing contribution to the local part 
of n comes from the third order diagram, which will give 
a constant. 

Now we go from the polarization to the lattice suscepti¬ 
bility ( 12 ) and expand it with respect to 1 - 

This is valid at high frequencies, since Xui decays as 
and we assume does not increase as a function of w, 
which will be justified a posteriori: 


^q,uj — Xi^ X<^B-q^cjX(^ iXc^ Xi^Bq^cjXcj) ^q)- 

(56) 


Applying the self-consistency condition (18 1 gives 

0 ^ ) Vq Xu} 

q 

0 =E [nq,c.-(l + X<.nq,^)'(A^-Fq) 
q 

q q 

-E(l + A-nq,^)"(A^-Vq). 


(57) 


According to the previous discussion, the first term van¬ 
ishes, the second term goes to a constant at large fre¬ 
quency, and 1 -I- also goes to a constant, so A^^ 

has to go to a constant as well. 


The nonzero retarded interaction at high frequencies is 
notably different from EDMFT and DB with the “dual” 
self-consistency condition. This frequency independent 
retarded interaction can be considered as a change in the 
instantaneous effective impurity interaction. Physically, 
this is the screening of the impurity by the collective 
modes. 


VIII. CONCLUSION 

In this paper, we formulated the inherently self- 
consistent dual boson scheme, capable of treating elec¬ 
tron structure and collective excitations in correlated sys¬ 
tems. The scheme employs the effective impurity prob¬ 
lem ([^ with fermionic and bosonic hybridization func¬ 
tions. We have shown that a proper choice of the retarded 
interaction leads to a significant simplification of the DB 
perturbation theory: all diagrams with local two-particle 
lines can be removed. In this way, the effective impurity 
problem contains more information about the nonlocal 
bosonic (e.g. charge, spin) fluctuations. Physically, the 
nonlocal bosonic fluctuations are determined both by the 
direct boson-boson interaction in the Hamiltonian and 
by the collective behavior of the fermions. These two 
phenomena need to be taken into account on the same 
footing. 

We have found that the nonlocal charge fluctuations 
make the system more insulating. However, to make a 
truly meaningful statement on this issue, the nonlocal 
spin fluctuations should be taken into account on the 
same footing, since they are expected to have a compen¬ 
sating effect. 

Perhaps somewhat surprisingly, an impurity model 
with a static U is not the optimal reference for the Hub¬ 
bard model with local interaction only. In finite dimen¬ 
sions, nonlocal charge and spin fluctuations can have 
large effectJ^^ and ideally the impurity problem should 
know about th is. DMFT overestimates the Neel temper- 
atur d^^ l ^^ l ^^ l ^^ l since the feedback of the strong antiferro¬ 
magnetic fluctuations on the impurity model is not taken 
into account. In 2d, these long-range antiferromagnetic 
fluctuations should bring down the Neel temperature all 
the way to zero according to the Mermin-Wagner theo¬ 
rem. 

In the presence of nearest-neighbor interaction V, we 
have found that the effective impurity interaction de¬ 
pends approximately linearly on V, as expected from the 
Peierls-Feynman-Bogoliubov variational principl^^. 

We have also presented a proof of the charge conserva¬ 
tion law in DB, and have shown that charge conservation 
can be achieved independently of the choice of retarded 
interaction. Numerically the DB charge susceptibility in¬ 
deed satisfies this law. Charge conservation is important 
for a proper description of long-wavelength modes. 

Finally, it is worth to draw an analogy to the DMFT 
paradigm. DMFT is known to be the best single¬ 
impurity method describing one-particle physics of cor- 
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related systems in a fully self-consistent way. It is im¬ 
portant that it uses only the single-particle Green’s func¬ 
tion of the impurity model. This is because of the proper 
choice of the fermionic hybridization function, which al¬ 
lows low-order dual fermion corrections to vanish (as 
these would otherwise include vertex parts of the im¬ 
purity). The theory presented above follows the same 
ideology for the collective excitations: the calculation re¬ 
quires knowledge of one- and two-particle properties of 
the self-consistent impurity model, whereas the proper 
choice of the fermionic and bosonic hybridization func¬ 
tions eliminates the diagrams with higher-order vertexes. 
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Appendix A: Sum rule for the fermion—boson vertex 


The fermion-bos on v ertex A and the fermion-fermion 
vertex 7 are relatecP*^, since they are both two-particle 
impurity correlators 



The impurity susceptibility is also a two-particle corre¬ 
lator, and it is related to the fermion-fermion vertex by 
a ladder equatioiP 

\uj ^ ^ 9v(y9v+uj^(7 ^ ^ QuaQv+u) loQv' tyQv'+UJ^(7• 

vv' go' 

(A2) 


Gombining these equations gives a useful sum rule for the 
fermion-boson vertex 

VG VV'gg' 

va 

= -x^ 

^ , 9iy<y9i'+i^ty^i/uj ~ 1 - (-^3) 

17(7 

This identity can also be obtained straightforwardly 
from the definition of A 

'y ' 9l7<79v+UJl7Xc^^l'UJ ~ y ^ {,^I7<7^U+U:<7'^^^'} 917(7 (^) 

17(7 17(7 

= (rit^nu,) - (n) (n) 

= -Xuj- (A4) 

Appendix B: Impurity vertex functions 

The three-particle vertex functions are given by 

772,3.4,5.6 = (5iS253545556)"^X 
[(C 1 C 2 C 3 C 4 C 5 C 6 ) - {cicl)-f‘^’°g2939596+ 

(ciCg) 7'^’°52ff35456 - 92939495 + 

(C 2 C 4 ) 91939595 - (C 2 C 5 ) 1 '^’°91939495+ 

{C 2 cl) 91939495 - (C 3 C 4 ) 91929595 + 

{C3cl) 91929495 - {C3cl) 91929495 - 

(C 1 C 4 ) (c 2 C 6 ) (C 3 C 5 ) -b (C 1 C 4 ) (C 2 C 5 ) {c^cD - 
{cicD (C 2 C 4 ) (C 3 C 6 ) -b (cic^) {C 2 cl) (C 3 C 4 ) - 
(ClCg) (C 2 C 5 ) (C 3 C 4 ) -b (ClCg) (C 2 C 4 ) (C 3 C 5 ) ] (Bl) 

whereas 7"*^’^ and 7^’^ can be expressed via 7®’° as follows: 

7i,2,5,6;3 = X3 2,3+4,4,5,653+4^4 + 

4 

c 4,0 c 4,0 I c 4,0 1 

472,34-4,5,6534-4 ~ ^3+4,574 2,4,654 + 534-4,67l,2,4,554j 

(B2) 


776;2,3 = +2X3) + 

E [7 1,’44-2,4,54-3,5,65455544-255+3 + 

4.5 

1^1,5^4+2.6 (C5+3C4) + 

1^1,4^5+3,6 (C4+2C5) + 
r 4.0 

01,475.4+2.5+3,65554+255+3- 
X 4.0 

01,574,6.4+2,5+3545654+255+3 — 

X 4.0 

01,475.6.4+2,5+3555654+255+3- 
1^6,5+37i,4.5.4+251545554+2 + 

1 ^ 6 , 4 + 271 , 4 . 5 . 5 + 351545555 + 3 ] • (B3) 
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This shows that 7 ^’^ contains one-particle and two- 
particle redncible contributions, whereas 7 ^’^ only con¬ 
tains the two-particle reducible part. 


Appendix C: Equivalence of the two self-consistency 
conditions within the EDMFT 


In EDMFT, where the dual polarization is equal to 


zero, Ilq^^ = 0 , the lattice susceptibility is given by ( 12 ) 

= + K-v,)-\ (Cl) 

Then, the “lattice” self-consistency condition is 

(C 2 ) 


At the same time, since Ilq = 0, the dual susceptibility 
X is equal to the bare dual susceptibility (|^. The old 
self-consistency condition (18) is 




(C3) 


and (C3) are clearly equivalent. This proves 


and 

the equivalence of the self-consistency conditions (|18|) 
and dlTl) in EDMFT. 


Appendix D: Relation between lattice and dual 
susceptibilities 


This expression can be illustrated by the diagrammatic 
series 

(D4)’ 


where the dashed line is the impurity susceptibility y, 
the filled ellipse denotes the interaction V — A and the 
triangles with fermion lines are the dual polarization 11 . 
Using the Eq. (D2| and Eq. (D3), one can get the re¬ 


lation between lattice and dual susceptibilities 


A = A 


X + xHy 

Xiy - K)x' 


(D5) 


Finally, we wish to get rid of U — A. To do this we solve 
Eq. (D2) for U — A and substitute the result 


V-A = 


X 


A(x + xnx) + x^ 

in Eq. (D5|. This results in Eq. ( [T^ 

- (x + xiix) [^(x + xiix) + x^] 


(D 6 ) 


A = A- 

xAy 

= (i + xn)[A(i + nx) + x]. 


(D7) 


The self-consistency condition 0 is formulated in 
terms of lattice quantities. We wish to have an interpre¬ 
tation of this self-consistency condition in the dual per¬ 
turbation theory, which requires us to rewrite the self- 
consistency condition in terms of dual quantities. To 
achieve this, we start by rewriting the bare dual bosonic 
propagator of Eq. 


[x-'-(^^-A)] '-X 
^ X(^-A)x 
1-(U-A)x' 


(Dl) 


The momentum and frequency labels in these expres¬ 
sion have been suppressed to simplify the notation. The 
Dyson equation gives the renormalized dual susceptibil¬ 
ity 


AO 

1 -AOfl 

xiy - A)x 

1 - (u-A)(x + xnx) 


(D2) 


The lattice susceptibility [Eq. 0 ] can be written in a 
similar way as 


X + xHx 

1- (u-A)(x + xnx) 


(D3) 


Appendix E: One— and two—particle reducible parts 


Here, we study the reducible parts of the higher-order 
vertices with respect to the purely fermionic vertex func¬ 
tions y'^’O and yO’O. We start with the simplest one- 
particle correction to the self energy, from the reducible 
part of Fig. where the one-particle reducible con¬ 
tribution comes from (B3). Just like 7 ^’^ = A, 7 ^’^ is 


nonzero even in noninteracting systems. Diagrammati- 
cally, the associated dual self-energy is 


-V®.. 

Sf 0 0 = 0^° 1'^ correction 

IX 




(El) 


where the dashed triangular vertex is the bare fermion- 
boson vertex 7 g’^ = — 1/x- This one-particle reducible 
part of acts as a correction to the dual self-energy in 
Eq.jSTJ 

The dual polarization for this case also has a simple 
form, 



This contribution to the polarization function corre¬ 
sponds to the limit U = A = 0. In that case, the in¬ 
teraction V is small and we can expand the relation for 
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the dual susceptibility to first order in V 

x(V-A)x 


X = 


1- (V - A)(x + x^x) 


= X^X- 


(E3) 


Since the local part of interaction Vq is equal to zero, 
the local part of the dual susceptibility is equal to zero 
as well. This means that the self-energy in (El) 

is equal to zero and therefore only bare fermionic lines 
remain in polarization function ( |E2[ ). 

In this case the lattice susceptibility i^ 


X = 


GG 


1-EGG 


(E4) 


This result is exactly the RPA relation for the susceptibil¬ 
ity. It means that for the simplest case the self-consistent 
DB approach reproduces the result of the Random Phase 
Approximation (RPA) method. In this limit of weak in¬ 
teraction, RPA indeed correctly describes the bosonic de¬ 
grees of freedonP^ and self-consistent dual boson reduces 
to RPA in the same limit. 

Now, we can go further and study the two-particle 
contribution from the reducible parts of the vertices 7 ^’^ 
and 7 ^’^, see Eqs. (B2) and (B3). Including the four- 
point vertex 7 "^’°, the additional reducible parts of the 
diagrams for the self-energy are 



yc _ 
^2PR — 




correction 


correction 


e 




(E5) 


When the impurity Green’s function is small, e.g., 
at high temperature and large U, these terms become 
small corrections. In the limit of small U, they also be¬ 
come small corrections since the three-particle correla¬ 
tors are small themselves. Thus, the main contribution 
of the three-particle correlation interaction functions to 
the self-energy comes from irreducible parts of 7 ^’^ and 
7 ^’^ and the vertex function 7 ®’° itself. 


Appendix F: Contribution of higher—order vertices 
to the self—energy 

A main advantage of the DB approach is the fact that it 
can be successfully applied to both the large and small U 
limits. For the small U case the diagrammatic expansion 
has a perturbative form and the higher-order vertices 
in U can be neglected. Here we analyze the physical 
irreducible contribution from the higher-order vertices, 
which are similar to those removed by the bosonic self- 
consistency condition. Other reducible contributions are 
also small for this case, as was shown in Appendix [E| 


The irreducible part of 7 ^’^ vertex function can be pre¬ 
sented in terms of the six-point vertex function as 

7i,’6;2,3 “ (A2X3) ^ 7 l.’4.4-r2.5.5-r3.6ff4g5g4-|-2g5-l-3- 

4.5 

(FI) 


One can rewrite this equation in following way 

7i. 6;2,3 X/ 54ff4+255ff5-l-3(l + ^ 74,4-|-2,4',4'+25'4'54'-|-2) X 
4,5 4' 

{^ + J2^^’°35'g5'+3) = XI'^74.4+2.5,5+3.65454-^25555+3, 
5' 4,5 

(F2) 


or equally 

X 7y6;2.3X2A35454+25555+3 X 
4,5 

( ~ 7 ~ X 74,4+2,4',4'+254'54'+2) X 

X2 X2 ^ 

/If 4,0 \ 

( T/ T/ 2^ 75,5+3,5',5'+355'55'+3j — 


X3 X3 

E 6,0 

71,4,4+2,5,5+3,65454+25555+3- 

4,5 


(F3) 


Using the relation (11) we get 

X 7y6;2.3X2X35454+25555+3X 
4.5 

/ 2 2,1 W 2 2,1 \ 

( -v 5'4,4+2;2j( 75,5+3;3j 

X2 X3 

E 6,0 

7l,4,4+2,5,5+3,65454+25555+3- 
4.5 


(F4) 


If we look again on the relation (11) we can realize that 
the bare vertex function — F can be expressed through 
the local three- and four-point vertex functions and local 
Green’s function: 


-1 

X 

= 7 "’" + 

“1 4 0 

— gg7 ’ 

X 

'1 

= < + 



(F5) 

The contribution of the six-point vertex to the self¬ 
energy E® now looks as follows: 


main term 


correction 



(F6) 


When the two-particle vertex or impurity Green func¬ 
tion g are small, then all terms with the dashed lines are 













16 


accounted for the small correction. Therefore we con¬ 
sider only the main term in the bare three-point ver¬ 
tex — ^ 74 4 _|_ 2.2 with the proper summation over the 

fermionic frequency. This is the case for example when 
the three-point vertex weakly depends on the fermionic 
frequency. Then one can obtain the following relation 


E 6,0 

'l'l,4,4-|-2,5,5+3,6ff45'4-|-23555-l-3 — 

4,5 

E 2 2 2 1 2 1 

'l'l’6;2, 3 X 2 X 374 ,’4+2;275’5-l-3;3ff454-l-2ff555+3- (F7) 

4,5 


Similarly to the six-point case we can obtain the relation 
for the five-point vertex as follows 


E 4.1 

7l,4,4-|-2,6;354ff4-|-2 — 
4 


7l’6;2,3X274,’4+2;354ff4+2, 


(F8) 


The expressions (F7) and (F81 are similar to those shown 
in Fig. and differ from them only by the replacement of 
the local Green’s functions by their dual counterparts. 
One can expect, however, that these diagrams, being 
physically important, indeed give the dominant contribu¬ 
tion to the self-energy, which therefore can be cancelled 
by the self-consistency condition, as it was shown in (291. 


Appendix G: Quadratic dependence on V in 
EDMFT 


Our EDMFT results show a quadratic dependence on 
V. This can be understood from a symmetry of the 
square lattice, and of hypercubic lattices in general. The 


Fourier transform of the nearest-neighbor interaction is 

d 

Vq = 2V COs{ka). (Gl) 

Q: = l 


All the momenta occur in pairs q, q', with q' = q-|- (tt, tt) 
and Vq = — Vqi. In the EDMFT formalism, Vq is the only 
quantity carrying momentum dependence, and it only 
occurs in formulas with a momentum-averaging, such as 
the self-consistency condition. Changing the sign of the 
interaction, i.e., Vq = —Vq does not change the self- 
consistency condition: 




Xw + — Vq 


= E- 


Xw + 


=E- 


= E-r 


X(^ + — [—Eq'] Xl^ + — Eq, 


(G2) 


This shows that the EDMFT result in a hypercubic lat¬ 
tice does not depend on the sign of the nearest-neighbor 
interaction. This is not what is expected physically, at¬ 
tractive and repulsive nearest-neighbor interactions are 
notably different. 

In self-consistent DB the results do depend on the sign 
of the interaction. Since there is an additional quantity 
with momentum dependence, Ilq^j, the sign flipping ar¬ 
gument of the previous paragraph no longer applies^ 

Diagrammatically, this can be seen from Eq. (D4|, 
where the momentum dependence comes in through the 
filled ellipse (via the interaction Eq) and through the 
fermionic bubbles (the dual polarization). EDMET does 
not contain the latter, so all of the momentum depen¬ 
dence comes from Eq. Due to the symmetry of the square 
lattice, changing the nearest-neighbour interaction from 
repulsive to attractive only rearranges the interaction 
terms and leaves their momentum averages the same. In 
dual boson, Eq. (D4) contains contributions where the 
both the dual polarization and the interaction appear. 
Exactly these diagrams give the contributions linear in 
E. 
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